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The strain in AlGaN/GaN heterojunction field-effect transistors (HFETs) is examined theoretically in the con- 
text of the fully-coupled equation of state for piezoelectric materials. Using a simple analytical model, it is 
shown that, in the absence of a two-dimensional electron gas (2DEG), the out-of-plane strain obtained without 
electromechanical coupling is in error by about 30% for an Al fraction of 0.3. This result has consequences for 
the calculation of quantities that depend directly on the strain tensor. These quantities include the eigenstates 
and electrostatic potential in AlGaN/GaN heterostructures. It is shown that for an HFET, the electromechanical 
coupling is screened by the 2DEG. Results for the electromechanical model, including the 2DEG, indicate that 
the standard (decoupled) strain model is a reasonable approximation for HFET calculataions. The analytical re- 
sults are supported by a self-consistent Schrodinger-Poisson calculation that includes the fully-coupled equation 
of state together with the charge-balance equation. 
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I. INTRODUCTION 

The piezoelectric and spontaneous polarization properties 
of the Al^Gai-a-N material system makes it attractive for cer- 
tain wide bandgap device applications. Heterojunction field- 
effect transistors (HFETs) made from AlGaN/GaN possess 
two-dimensional electron gas (2DEG) concentrations on the 
order of 10 13 cm~ 2 at the AlGaN/GaN interface. Owing to 
the polarization effects, the 2DEG can be established at zero 
gate bias and achieved without intentional doping. 1-6 A back 
doping design 7 can further increase the 2DEG concentration 
to approximately 3 x 10 13 cm -2 . The dense 2DEGs result in 
r.f. device power densities in the range 6-12 W/mm, making 
AlGaN/GaN HFETs attractive for high-power applications. 

The formation of the 2DEG can be traced directly to the 
large built-in polarization fields in the material system. In 
addition to the piezoelectric polarization due to the strain in- 
duced by the lattice mismatch, there is a spontaneous polar- 
ization in the direction N — > Ga along the c-axis which can be 
several times larger than the piezoelectric polarization. These 
polarizations induce local electric fields which, in turn, cause 
a notch to form in the conduction band edge at the interface. 
Electrons congregate in the notch, either from surface states 
or from unintentional dopants and traps, forming the 2DEG. 

In previous HFET analyses, the electrical and mechanical 
properties of the epitaxial layers have been treated indepen- 
dently. The standard elastic theory is applied and Hooke's 
law is assumed to hold in obtaining the strain tensors in each 
material layer. In effect, this assumption decouples the strain 
tensor from the Poisson and Schrodinger equations and ren- 
ders the task of computing the free charge distribution and 
band structure more tractable. In reality, however, the elec- 



trical and mechanical properties of piezoelectric materials are 
coupled as seen from thermodynamics 8 and are best treated 
using a unified approach, especially since the piezoelectric re- 
sponse is large in AlGaN for high Al fractions. 

In other work, large modifications to the elastic prop- 
erties have been predicted for buried group Ill-nitride 
nanostructures 9 and free-standing plates of A1N. 10 The lat- 
ter is a geometry closely related to the AlGaN barrier layer 
of the HFET. As a direct comparison, we take an A1N plate 
experiencing biaxial strain equivalent to the strain of a pseu- 
domorphic A1N layer on GaN, or r y X x=lyy=— 0.024. The re- 
sulting strain using an uncoupled calculation in the absence of 
spontaneous polarization is 7 22 =0.014 and the electric field 
is E z =— 5.1 MV/cm. In the coupled calculation, we obtain 
7^=0.012 and E z =— 4.7 MV/cm. Based on these substan- 
tial differences there is strong motivation to make a detailed 
analysis of a transistor structure including a detailed analysis 
of the spontaneous polarization moment and free carrier dis- 
tributions. 

In the present work, the strain and electric fields obtained 
from the fully-coupled model are contrasted with those from 
the standard model for a representative HFET structure. The 
deviations from the standard model are discussed in detail. 
The theory uses the fully-coupled piezoelectric equation of 
state to obtain the strain field. A simple analytical model is 
then used to obtain the interaction between the polarization- 
induced electric field and the strain tensor. Also included is 
the effect of 2DEG charge-screening on the electromechani- 
cal coupling. The results from the analytical model are sup- 
ported by those from a self-consistent Schrodinger-Poisson 
model that includes charge balance and free hole states. 

This paper is organized as follows: In Sec. II, the strain 
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tensor is derived from the fully-coupled equation of state. 
The Poisson equation is solved analytically for a bi-layer Al- 
GaN/GaN slab and the strain components and electric field 
in the barrier material is derived in closed form. In Sec. Ill, 
the calculated deviations of the standard model relative to the 
fully-coupled model are shown. The effect of screening from 
the 2DEG is discussed. The results from a fully-coupled nu- 
merical model are presented. The results are summarized in 
Sec. IV. 



n. MODEL DESCRIPTION 

The fully-coupled equation of state for piezoelectric mate- 
rials is expressed as 11 



(1) 



in component form where is the stress tensor, Cijki is the 
fourth-ranked elastic stiffness tensor, jki is the strain tensor, 
ekij is the third-ranked piezoelectric coefficient tensor, is 
the electric field vector, and the indices i, j, k, and I run over 
the Cartesian coordinates x, y, and z. Summation over re- 
peated indices has been suppressed. Eq. (1) in the absence 
of piezoelectric coupling (ekij— *-0) is recognized as Hooke's 
law relating stress to strain. The fully-coupled model is now 
applied to the representative HFET structure shown in Fig. 1. 
A one-dimensional (1-D) approach is used where the strain 
in the bi-layer slab is assumed to be piecewise homogeneous. 
This approximation means that 7^ is diagonal for the coordi- 
nates defined in Fig. 1 and that j xx and -f yy are, to first order, 
uniform in each layer. The strain induces a piezoelectric po- 
larization given by 



P? 



(2) 



From straightforward electrostatics we know that Eqs. (1) and 
(2) are coupled through the Poisson equation which is written 
in our notation as 



V-eV<£= -p + V- (P s + P p ), 



(3) 



where <j> is the electrostatic potential, e is the electric permit- 
tivity, p is the charge distribution and includes both ionic and 
free charges, P' s is the spontaneous polarization, and P p has 
been defined in Eq. (2). In the absence of shear strains, P s and 
P p are directed in the [0001] direction for a cation-faced struc- 
ture as shown in Fig. 1, and with the sixfold rotational symme- 
try along the c-axis, 7x X =7 OT/ . From Eq. (2), the piezoelectric 
polarization along [0001] for wurtzite materials then becomes 



P p = 2e 3 i7 xx + e 337z 



(4) 



It is assumed that j xx is determined by the atomic alignment 
in the c-plane across the AlGaN/GaN interface. This pseudo- 
morphic boundary condition gives 



a b 



in the AlGaN layer and "/ xx —0 in the GaN layer, where a a 
and ai, are the c-plane lattice constants in the AlGaN and GaN 
layers, respectively. In reality, besides the force exerted on 
the AlGaN layer by the GaN layer, one would also expect the 
AlGaN layer to exert a force on the GaN layer, causing some 
strain on the GaN side of the interface. This local distortion of 
the lattice would be a maximum at the junction and become 
vanishingly small towards the interior of the GaN. Such re- 
finements, however, would require a three-dimensional (3-D) 
modeling of the electromechanical problem, as it would have 
to include non-zero and position-dependent shear terms in the 
strain tensor. This is beyond the scope of the present work and 
is not essential for our investigation. 

Thus, within the present work, ^ xx is assumed to be de- 
coupled from the electrical properties, essentially fixed by the 
bonding arrangement at the interfacial plane, with "f zz serving 
as the principal vehicle for electromechanical coupling. From 
Eq. (1) and from the foregoing assumptions, the normal com- 
ponent of the stress tensor is given by 



a zz = 2Ci37 xx + Cw-fzz - e 33 E z 



(6) 



From the boundary conditions for a free surface, the stresses 
along the outward normal must be zero. As a consequence, 
Eq. (6) then gives for the 1-D case 
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which is the out-of-plane strain, where E z —~dc/)/dz. From 
Eq. (7), it is seen that the first term is consistent with the Pois- 
son ratio of the standard model for j zz and a new coupled term 
representing the electromechanical coupling. To date, the lat- 
ter has not been included in AlGaN/GaN HFET models. The 
extent of the coupling will depend on the Al mole fraction and 
the electric field in the layers, which, in turn, will depend on 
the charge distribution and boundary conditions used to solve 
Eq. (3). The actual problem is quite complicated and involves 
the self-consistent solution of the equation of state in conjunc- 
tion with the Poisson, Schrodinger, and charge balance equa- 
tions. To illustrate the concept of electromechanical coupling 
without obscuring much of the underlying physics, a simple 
example has been contrived from which analytical results can 
be extracted so that the effect of coupling can be easily quan- 
tified. The charges assumed herein are the polarization space 
charges and a sheet electron distribution at the AlGaN/GaN 
interface. Ordinarily, the electron distribution should be cal- 
culated quantum-mechanically from the electron eigenstates, 
as described in Ref. 12. In order to illustrate these effects an- 
alytically, the 2DEG is modeled here as a <5 function localized 
at the AlGaN/GaN interface. These induced electrons may 
be accounted for by surface charges such as donor states, by 
deep-level traps in the AlGaN and GaN layers, or by unin- 
tentional dopants. The exact origin of the induced 2DEG is 
still unknown, but its presence is well-established from pre- 
vious work. 1 As will be seen later, representative values for 
the induced 2DEG are used. The 1-D Poisson equation then 
becomes 
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where is the areal 2DEG concentration and e Q is the mag- 
nitude of the electronic charge. Equation (8) is integrated to 
obtain subject to the boundary conditions 0=0 at z~0 and 
z~t a +tb and also to the continuity of and the electric dis- 
placement across the interface. The condition = at the 
two surfaces presumes that the polarization charge of the ma- 
terials is terminated by external charges and that the applied 
bias voltage is zero. The general solution in each layer, ob- 
tained by integrating Eq. (8) twice, is given by 

, ^ s , 2(e 31 C 3 3 ~ e 33 C 13 ) lxx A 

4>=—z+ — z + —z + B, (9) 

e e'C 33 e' 



where 



-33 



(10) 
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and is layer-dependent because of the varying elastic and 
piezoelectric coefficients, and A and B are unknown con- 
stants. Thus there are four unknowns, two in each layer. The 
B's are eliminated by enforcing the boundary conditions 0=0 
at z=0 and z=t a +tb. It should also be noted that -f xx = in 
the GaN buffer because of the simplifying assumptions made 
previously. The relationship between the two ^4's is estab- 
lished from the continuity of the electric displacement and is 
found by integrating Eq. (8) across the AlGaN/GaN interface 
(see Fig. 1) encompassing the 2DEG. This boundary condi- 
tion gives the relation between the electric fields in the two 
layers as 



dz 



'° _ p S |*i" + 2 ( e 3lC 3 3 - 633613)7^ 
f~ ta C33 



+ e Q n 
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(11) 

From Eq. (1 1), we obtain the relation between the A's in the 
respective layers as 



^GaN _ ^AlGaN 



,(2) 



(12) 



From Eq. (12) and from the continuity of across the inter- 
face, we can solve for all of the remaining unknowns. 

The results are greatly simplified in the limit tt,^>t a , a con- 
dition that is met for most HFETs. In this limit, the electro- 
static potential in the AlGaN layer is given by 



CAlGaN 



ps 
AlGaN 



GaN , 
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-Z + 



e 3lC33 — e 33 Ci3 



(13) 



and is constant in the GaN layer. From Eq. (13), the electric 
field in the AlGaN layer is given by 



The electric field within the standard model is obtained by 
replacing e' with e in Eq. (14). It is seen that for growth 
along [0001], the electric field in the AlGaN from the coupled 
model is smaller that of the standard model, since e 2 3 /C33 
is always positive. The vertical strain component within the 
fully-coupled model is obtained from 



Izz = 



-2— 7^ + 
(^33 



2e33(e33Ci3 — 631(733) 



Ixx + 



C 33 (eC 33 + el 3 ) 

e33(P<L N -^AlGaN + ^ (2) ) 



eC 33 + el 



(15) 
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The first term in Eq. (15) is the typical result obtained from 
uncoupled elastic theory. The remaining terms are due to elec- 
tromechanical coupling. 

Besides the simple analytical approach described in this 
work, a previously-published numerical Schrodinger-Poisson 
model 12,13 has been modified to include full coupling. The 
elements of the strain tensor enter the Schrodinger equation 
via deformation potential theory. One of these elements, -f zz , 
is now a function of the electrostatic potential calculated from 
the Poisson equation. In turn, the electrostatic potential is a 
function of the piezoelectric polarization, itself a function of 
Izz- With each Schrodinger-Poisson iteration, the strain must 
be updated via Eq. (7). Also included is the charge-balance 
equation as described in Ref. 12. The calculation is regarded 
as having converged when the maximum change in the elec- 
trostatic potential between the current and the previous itera- 
tion is less than 0.01 meV. 



III. RESULTS AND DISCUSSION 

This section contains calculations used to compare the un- 
coupled and coupled models for a model HFET structure. The 
presence of the 2DEG and its effect on the problem are investi- 
gated. Table I shows the material parameters 14-20 used in the 
calculations. The signs of the polarization parameters are de- 
fined in relation to the [0001] direction: a negative sign means 
that the vector is in the [0001] direction. 

The calculated strain for both the standard and fully- 
coupled models as a function of the Al fraction of the barrier 
layer is shown in Fig. 2. It is seen that the electromechanical 
coupling changes j zz significantly in the absence of a 2DEG. 
It is also noted that the coupling reduces ^ zz relative to the 
standard model. This result is evident from Eq. (7) in which 
it is seen that the coupling opposes the contraction of material 
along the c axis when there is tensile strain in the c plane. A re- 
alistic comparison between the two approaches should include 
the mobile charge drawn to the interface by the polarization 
space charge. The net polarization space charge at the inter- 
face is given by P=[Pg^- P mg^- p I\g^I and for a 30% 



Al alloy composition is approximately 1.63 x 10 13 e Q /cm 2 . 
The 2DEG can be expected to neutralize the space charge to 
some extent. In Fig. 2, is taken to be 0.8P/e o , although it 
will be seen shortly from a Schrodinger-Poisson solution that 
the neutrality factor is closer to 0.9. With screening from the 
2DEG included, the discrepancy between the two models is 
reduced. This result is due to the reduction in the polarization- 
induced electric field caused by screening and is shown more 
clearly in Fig. 3. The magnitude of the 2DEG chosen here is 
representative of the induced electrons seen in actual HFET 
devices. 1 We do not, however, explore the origin of the 2DEG 
in this simple approach, as the only purpose in including the 
2DEG is to show that electron screening has a strong influence 
on the electromechanical coupling. 

Next we examine the deviation A between the coupled and 
uncoupled models. The deviation in the out-of-plane strain 
component ~f zz is defined as 

uncoupled 

A = lzz - lzz (16) 

where ^ zz is given by Eq. (15) for both the screened and un- 
screened cases. The calculated A is shown in Fig. 4. For 
a typical HFET with an Al mole fraction of 0.3, the error is 
about 30% if the effect of the 2DEG is neglected, as compared 
to 17% for free-standing plates of A1N. It is clear that if the Al- 
GaN layer is considered in isolation, the standard model gives 
a large error in the strain tensor. In an HFET structure, how- 
ever, the importance of the AlGaN/GaN interface cannot be 
overestimated. This interface introduces additional physics, 
the most important being the 2DEG. Effectively, the presence 
of a large 2DEG restores the validity of the standard model. 

Fig. 5 shows the difference in j zz when the spontaneous po- 
larization is excluded. It is well-known that the spontaneous 
polarization is the larger of the two polarizations in the Al- 
GaN/GaN material system. For example, for a Al fraction 
of 0.3, the space charge induced by the discontinuity of the 
spontaneous polarization is about 9.8 x 10 12 e Q /cm 2 as com- 
pared to about 6.7 x 10 12 e Q /cm 2 for the piezoelectric po- 
larization. The difference increases for higher Al fractions. 
Accordingly, with just the piezoelectric polarization included, 
the electromechanical coupling is underestimated. This is be- 
cause the electric field is reduced compared to Fig. 3. Once 
again, however, the screening from the 2DEG reduces the cou- 
pling. 

Fig. 6(a) shows the conduction band edge and the free 
electron distribution obtained from the Schrodinger-Poisson 
calculation that includes the full electromechanical coupling. 
The structure is a standard HFET design consisting of 300 
A of Alo.3Gao.7N on a thick GaN buffer. The 2DEG is as- 
sumed to come from surface donor states 12 having an acti- 
vation energy of 1.4 eV below the conduction band edge at 
the surface, effectively pinning the surface potential at 1 .4 V. 
The resulting 2DEG concentration after convergence is about 
1.47 x 10 13 cm~ 2 . This magnitude is about 0.9P/e o , almost 
neutralizing the polarization-induced space charge. As an ex- 
ample, the electric field in the Alo.3Gao.7N layer in the ab- 
sence of a 2DEG is about 3.1 MV/cm from Fig. 3. The fully- 



coupled Schrodinger-Poisson model gives a field of about 0.45 
MV/cm and the standard Schrodinger-Poisson model a field of 
about 0.46 MV/cm. Within numerical errors, the 2DEG is al- 
most unchanged between the coupled and standard models as 
shown more clearly in Fig. 6(b). Consequently, there is little 
change in the conduction band edge between the two mod- 
els. The calculated eigenstates from the Schrodinger equation 
and the electrostatic potential from the Poisson equation show 
little change relative to the standard model after both models 
have converged to a solution. 



IV. SUMMARY AND CONCLUSIONS 

In summary, a fully-coupled electromechanical model has 
been presented for the strain in AlGaN/GaN HFETs. The 
model goes beyond the standard Hooke's law description and 
uses the fully-coupled equation of state for the stress-strain re- 
lationship. Using a simple analytical model for the HFET, it is 
shown that the present model gives rise to significant changes 
in the strain field and, consequently, to the electronic prop- 
erties of the device in relation to the standard model in the 
absence of free charges. When free charges are present in the 
form of a 2DEG, the built-in electric field is reduced through 
the screening of the fixed space charge. This effect, in turn, 
reduces the electromechanical coupling and brings the results 
from the standard and fully-coupled models into closer agree- 
ment. 
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TABLE I. Strain-related material parameters used in the present 
model. The elastic stiffness constants are in units of GPa and the 
piezoelectric stress constants and spontaneous polarization in units 
of C/m 2 . 

Material C\ 3 C33 a (A) e 3 i e 3 3 -Pspom e/ep 
ATN 108 a 373 a 3.1 12 b -0.58 d 1.55 d -0.081* 8.5 g 
GaN 103 a 405 a 3.189 c -0.36 c 1° -0.029* 10 s 



a Reference 14. 
b Reference 15. 
c Reference 16. 
d Reference 17. 
"Reference 18. 
f Reference 19. 
s Reference 20. 

h Contracted index notation: Ci3—C xxzz =C yyzz —C zzxx =C : 
Cs3=C zzzz , e 3 i—e zxx =e zy y, and ezz=e zzz . 
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FIG. 1. Cross-section of a model cation-faced AlGaN/GaN HFET 
showing the direction of the piezoelectric and spontaneous polariza- 
tion vectors in relation to the z-axis. t a and tb are the thicknesses of 
the AlGaN and GaN layers, respectively. 
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FIG. 4. Deviation of the out-of-plane strain for two coupled cases, 
with and without free-carrier screening. 
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FIG. 2. Calculated out-of-plane strain for the standard and cou- 
pled models. Two coupled cases are shown, one without free elec- 
trons and one with free electrons. 
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FIG. 5. Deviation of the out-of-plane strain for two coupled cases, 
with and without free-carrier screening with the spontaneous polar- 
izaion set to zero. 
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FIG. 3. Electric field in the AlGaN layer for the standard and cou- 
pled models with and without 2DEG screening. For the screened 

cases, n y 



- {2) = 0.8P. 
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FIG. 6. (a) Calculated conduction band edge (left y-axis) and free 
electron distribution (right j/-axis) along a c-axis slice of a repre- 
sentative HFET structure at 300 K. The Alo.3Gao.7N cap layer is 
300 A thick. The dot-dashed line is the Fermi energy. The surface 
of the structure is at the left extremity and the interior towards the 
right. The the 2DEG is assumed to originate from surface states. 
The standard (solid line) and coupled (dashed line) electron distribu- 
tions are shown, (b) The free electron distribution in (a) plotted over 
a narrower region of the structure to show the difference between the 
coupled and standard cases more clearly. 
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